Kolmogorov-Arnold-Moser— Renormalization-Group Analysis of Stability in 

Hamiltonian Flows 

M. Govin, C. Chandre, and H. R. Jauslin 

Laboratoire de Physique, CNRS, Universite de Bourgogne, BP 400, F-21011 Dijon, France 



00 

in 

(N 



(N . 
O ■ 
(N ■ 

o : 

00 , 

o\ ■ 
>v 

■ 

i ' 

o : 

^ . 

O ■ 



We study the stability and breakup of invariant tori 
in Hamiltonian flows using a combination of Kolmogorov- 
Arnold-Moser (KAM) theory and renormalization-group tech- 
niques. We implement the scheme numerically for a family of 
Hamiltonians quadratic in the actions to analyze the strong 
coupling regime. We show that the KAM iteration converges 
up to the critical coupling at which the torus breaks up. 
Adding a renormalization consisting of a rescaling of phase 
space and a shift of resonances allows us to determine the 
critical coupling with higher accuracy. We determine a non- 
trivial fixed point and its universality properties. 

PACS numbers: 05.45. +b, 02.40.-k, 64.60.Ak 



Regular Hamiltonian dynamics is characterized by the 
existence of as many independent conserved quantities as 
degrees of freedom d. As a consequence, each trajectory 
is confined to evolve on an invariant torus of dimension d. 
This global regularity is destroyed in general by any small 
perturbation: some of the invariant tori break up and 
chaotic trajectories appear. However, many invariant tori 
survive and are just slightly deformed. The persistence 
of these tori has strong implications on the global phase 
space structure of the dynamics. They strongly limit the 
spreading of chaos and prevent global ergodic behavior. 
In particular, for two degrees of freedom, they block com- 
pletely the long distance diffusion in phase space. In this 
context, it is important to understand how invariant tori 
are destroyed as the size of the perturbation is increased. 

The Kolmogorov-Arnold-Moser (KAM) technique was 
developed originally to prove the stability under small 
perturbations of a large fraction of the invariant tori of 
an integrable Hamiltonian. It has been later applied to a 
very large spectrum of phenomena both in classical and 
in quantum mechanics. The KAM theorem (l]-D states 
that, if the frequency satisfies a diophantine condition, 
and the size of the perturbation e is smaller than an 
Eo ; then a torus of that frequency will be stable. The 
proof is based on an iterative algorithm to construct the 
invariant tori. Each step of the KAM iteration consists 
of a coordinate transformation that reduces the size of 
the perturbation from order e to e 2 . 

However, it was remarked early that the value £q ob- 
tained in the proofs (t]], for which the iteration is shown 
to converge, is much smaller than the critical coupling e c 
at which the torus is known 0] to become unstable and 
to break up into a Cantor set (Aubry-Mather set) 
Since there is no physical phenomenon associated with 



£o, the difficulty to prove the convergence of the KAM 
iteration for e e [eo,£ c ) was interpreted as a shortcom- 
ing of the iterative approach. A first goal of the present 
Letter is to provide numerical evidence that the KAM 
iteration can be expected to converge all the way to the 
critical coupling e c . 

Renormalization-group (RG) ideas were proposed for 
the analysis of the breakup of invariant tori ]| jll]] 
based on the observation of self-similar scaling proper- 
tics [^2|-p^[|, and inspired on the renormalization theory 
of phase transitions JL5| |. For Hamiltonian flows, Es- 
cande and Doveil set up an approximate renormalization 
scheme that consisted of a combination of KAM 
transformations with a rescaling and shift of resonances. 
The approximation consisted in retaining at each step 
only two Fourier modes, corresponding to the main reso- 
nances, and neglecting all the other modes as well as all 
terms of order higher than two in the momentum. Al- 
though the approximation seems very drastic, since the 
effect of the neglected terms could not be controlled, the 
general ideas of the renormalization group approaches 
suggest that if the essential features are kept in the ap- 
proximate scheme, the neglected parts may turn out to 
be irrelevant. This can of course only be firmly justified 
by treating the problem without the approximation, and 
establishing which degrees of freedom are relevant and 
which ones are irrelevant. 

An outstanding question in this context has been 
whether this approximate renormalization scheme can be 
improved in a systematic way, in order to understand 
the properties of an exact renormalization scheme. The 
idea of the renormalization approach is to iterate a trans- 
formation in the space of Hamiltonians. For couplings 
e < £ c the iteration should converge to a trivial fixed 
point. This can be considered as an alternative version 



of the KAM theorem |16|. The set of Hamiltonians with 
critical coupling e c form a surface in the space of Hamil- 
tonians. The renormalization transformation maps this 
surface into itself. The main hypothesis of the renormal- 
ization group approach is that there should be another 
nontrivial fixed point on the critical surface that is at- 
tractive for Hamiltonians on that surface. Its existence 
has strong implications concerning universal properties 
in the mechanism of the breakup of invariant tori. 

The analysis of the approximate scheme of Escande 
and Doveil HD, and conceptually similar studies on the 
dynamics of area-preserving maps Ej, give support to 
the validity of the general picture. In order to give a 



1 



Phys. Rev. Lett. 79, 3881 (1997). 



2 



firm basis to this approach, some essential points need to 
be established: since the renormalization transformation 
contains at each step a KAM transformation, it is indis- 
pensable to rely on the existence of a well defined KAM 
transformation for e all the way to e c . Furthermore the 
KAM iteration should either converge for all e < e c or at 
least it should do so when combined with the rescaling 
and resonance shift. 

Our approach has been to implement the KAM algo- 
rithm numerically for Hamiltonians with two degrees of 
freedom, quadratic in the action variables A = (Ax,A2), 
of the form 

H(A, <p) = \m{<p){to ■ A) 2 + [w„ + eg(<p)n] ■ A + ef(<p) , 

(1) 

where ip — (ipi, ^2) are two angles, wo is the frequency 
vector of the considered torus, and fi is some other con- 
stant vector. This type of model can describe, e.g., a 
charged particle interacting with two waves in a plasma 
and has a wide range of applications. 

First we study the convergence properties of the KAM 
iteration. Then we define a renormalization transforma- 
tion (KAM-RG) by combining it with a shift of reso- 
nances, and a rescaling of momentum and energy. 

We show that the KAM algorithm can be used to deter- 
mine the critical coupling e c . We then show that this can 
also be achieved with the KAM-RG transformation, and 
that it is much more efficient, in that very high precision 
can be obtained already with few Fourier coefficients. 

We establish the existence of a nontrivial fixed point 
by iterating the KAM-RG transformation on the criti- 
cal surface. We characterize the nontrivial fixed point 
Hamiltonian by its Fourier representation with a very 
large number of coefficients (around 4000). 

We calculate the critical exponents by linearizing the 
renormalization map around the fixed point. The result 
shows that there is only one unstable direction, trans- 
verse to the critical surface. 

KAM iteration. — We consider a Hamiltonian of the 
form (Q), where eg(tp)tt and ef(if) are the perturbations. 
The KAM procedure consists of eliminating the pertur- 
bation of order e by a canonical transformation, which 
in turn produces terms of order e 2 . In the approach we 
take, following W. Thirring |L7|,[l8), m((p) is of order one 
and does not need to be eliminated. The point is that 
we want to establish the existence of a torus of frequency 
u>o- If (0) can be transformed into a Hamiltonian of the 
same form but with / = 0, g = 0, with arbitrary m(ip), 
the equations of motion show immediately that at A = 
there is an invariant torus of frequency u>o- The great 
advantage of this approach is that the KAM transforma- 
tion can be taken such that the transformed Hamiltonian 
is also quadratic in the actions, i.e. the KAM transfor- 
mation is a mapping among Hamiltonians of the form (p]) 



(f2 is not changed). This is very convenient for a numer- 
ical implementation, since at each step of the iteration 
the Hamiltonian is completely determined by three func- 
tions /, g, m of the angles ip, that can be represented by 
their Fourier coefficients. Thus the only approximation 
involved is the representation of the functions by a finite 
number of coefficients. In our implementation we take 
all the Fourier modes with wavenumber k € Z 2 in the 
square Cl of length 2L + 1 centered at (0, 0). 

We perform the KAM transformation as a canonical 
change of coordinates (A,(p) 1— > (A',cp') in the Lie rep- 
resentation p9| : 

H'(A',<p')=e+^ A ^H(A,cp) \ (AW) 

= H + e{S,H}+ £ -{S,{S,H}}... (2) 

generated by a function 5 1 , linear in the actions A, of the 
form 

S(A,(p) = Y{ip)Sl ■ A + Z{if) + a£l ■ <p . (3) 

The functions Y((p), Z(cp) and the constant a are chosen 
such that the terms of order e cancel out. We start with 
the same initial Hamiltonian as in Refs. 0(9): 

H{A,v) = \{Sl-A) 2 +u, -A + ef{v) , (4) 

where ft = (1,0), u> = (1/7,-1), 7 = (1 + V5)/2, 
and f(<p) — cos(i>'i • <p) + cos^ ■ if), where V\ — (1,0) 
and V2 — (1, 1). We perform an iteration of the above 
KAM transformation, representing all the functions by 
their Fourier series truncated by retaining only the co- 
efficients in the square Cl- For fixed L we take succes- 
sively larger couplings e and determine whether the it- 
eration converges to a Hamiltonian with / = 0, g = 
or diverges. By a bisection procedure we determine the 
breakup coupling e c (L). We repeat the calculation with 
larger numbers of Fourier coefficients, to obtain a more 
accurate approximation. In Fig. 1 we show s c (L) i.e. 
the dependence of the breakup coupling on the number 
of Fourier coefficients retained. We observe that e c (L) 
decreases with L in a stepwise manner. It stays essen- 
tially constant except at the points where a new rational 
approximant of the frequency loq is included, correspond- 
ing to a resonance at the next smaller scale. The size of 
the jumps diminishes approximately geometrically, and 
we can extrapolate to obtain the value e c (L) — * 0.0276. 
This value is close to the breakup e c = 0.0275856 ob- 
tained by the Greene criterion ^,^] which is surmised 
to yield the exact value. 

The calculation provides an explanation of the fact 
that this global convergence is not easy to establish an- 
alytically. The KAM theorem is proven by showing that 
for e < £0 the iteration is contractive. We observe in 
the numerical implementation of the iteration that for 
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larger e the convergence is not monotonic: in the first 
KAM steps the norm of the perturbation can actually 
grow and oscillate, until it reaches a region where it is 
contractive. 

KAM-renormalization group transformation. — In 
this section we show that by combining the KAM trans- 
formation with a rescaling and a shift of resonances we 
can greatly improve the determination of the breakup 
coupling. The intuition is that the renormalization treats 
specifically the effect of the main resonances at each 
scale. Their Fourier modes correspond to the rational 
approximants of u>o, which produce the small denomina- 
tors in the perturbation expansion and are responsible 
of the eventual breakup of the torus. Our results show 
that renormalization flattens out the step structure of 
the jumps in e c (L), and gives a very accurate value for 
the breakup coupling, already for a representation with 
a small number of Fourier coefficients. 

The KAM-RG transformation is composed of four 
steps: (1) a KAM transformation as described above, (2) 
a shift of the resonances: a canonical change of coordi- 
nates that maps the next pair of daughter resonances of 
the sequence of rational approximants into the two main 
resonances, (3) a rescaling of energy (or equivalently of 
time), and (4) a rescaling of the action variables (which 
is a generalized canonical transformation). 

The aim of this transformation is to go from one scale 
to a smaller scale. The last three steps are implemented 
as follows: The two main resonances (1, 0) and (1, 1) are 
replaced by the next pair of daughter resonances (2, 1) 
and (3, 2). This change is done via a canonical transfor- 
mation (A, ip) i-> (N~ 2 A, N 2 ip) with 



N 2 = 



2 1 
1 1 



This linear transformation multiplies ujq by 7 -2 ; there- 
fore we rescale the energy by a factor 7 2 in order to keep 
the frequency at o)o . Then we perform a rescaling of the 
action variables: we change the Hamiltonian H into 



H'(A,<p) = \H[j,<p 



with A such that the mean value of m remains equal to 1. 
This magnifies the size of the daughter resonances, and 
places them at the location of the original ones. 
The functions m, g, f are rescaled and ft = (l,a) is 
transformed as 



1 



a 



a 



2 + a 



(5) 



An effect of the shifting and rescaling is thus to select 
a = 7 _1 : Jl converges under successive iterations to O*, 
which is orthogonal to u>o and is the unstable eigenvector 
of N 2 with the largest eigenvalue j 2 . 



In Fig. 1 we show the values of the breakup coupling 
e c (L), calculated by this KAM-RG transformation. We 
obtain e c G [0.02758,0.02760], which is in very good 
agreement with the value e c = 0.0275856 obtained with 
the Greene criterion. 

The improvement with respect to the KAM iteration 
is not only quantitative; the disappearance of the steps is 
a strong evidence that the KAM-RG transformation we 
have constructed captures the essential physical mecha- 
nism of the breakup of the tori. 

Nontrivial fixed point. — By iterating the KAM-RG 
transformation starting from a point on the breakup sur- 
face, we observe that the process converges to a non- 
trivial fixed point ii*, which we characterize by the 
Fourier coefficients of the three functions f*, g*, and 
f2* = (1,7 _1 ). Figure 2 shows the weight of the Fourier 
coefficients of /*. We observe that the nonzero coeffi- 
cients are strongly concentrated on a band around the 
direction D± of the frequency vector ljq, i.e., perpendic- 
ular to the line of resonances. The decrease of the size 
of the coefficients along D± is quite slow. The Fourier 
coefficients of g* and m* have a similar overall behavior, 
but they decay faster in the D± direction. By linearizing 
the KAM-RG transformation around if*, we calculate 
the critical exponents. There is only one with \S\ > 1. 
This implies that the critical surface, which is the stable 
manifold of , is of codimension one. The value we ob- 
tain for the relevant critical exponent is 6 € [2.67,2.68] 
which is quite close to the one obtained by MacKay for 



area-preserving maps 11 (8 = 2.65), and to the one ob- 
tained by Escande et al. with the approximate scheme 
(5 = 2.75) @. 

We also obtain an attractive cycle of period three, as it 
had also been encountered in area-preserving maps p2[ . 
If odd perturbations are included, higher period can ap- 
pear. In the Hamiltonian case these cycles are simply 
related to the fixed point £/* by symmetry p4| . 

In conclusion, we have shown that the KAM-RG tech- 
nique can be implemented numerically with high accu- 
racy. The results show that the KAM-RG transforma- 
tion follows closely the mechanism of the breakup of the 
invariant tori by a sequence of resonances. Their effect 
at all scales leads to the universal self-similar structure 
of the critical tori. We have presented results for two de- 
grees of freedom, but the extension to three p3] or higher 
dimensional systems should be accessible. The extension 
of the scheme to other frequencies that are quadratic 
irrationals is relatively clear, but the case of a general 
irrational frequency will involve some qualitatively new 
features. 
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letti for providing the program to determine the critical 
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FIG. 1. Critical coupling e c {L) as a function of L, the 
size of the cell Cl containing (2L + l) 2 Fourier coefficients. 
The upper curve corresponds to the KAM transformation, 
and the lower one (enlarged in the inset) to the KAM-RG 
transformation. 




FIG. 2. Weight of the Fourier coefficients of /*: 
White: < 10" 10 , grey levels: [1(T 10 , 1CT 7 ], [1CT 7 , 
[lCT 5 ,lCr 3 ], black: [icr 3 ,i(r 2 ]. 
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